One-Dimensional Potential Model for Image States on Free-Standing Graphene. 
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In the framework of the non-local dielectric theory the static non-local self-energy of an electron 
near an ultra-thin polarizable layer has been calculated and applied to study image-states near free- 
standing graphene. The corresponding series of eigenvalues and eigenfunctions have been obtained 
by solving numerically the one-dimensional Schrodinger equation. We compare with the Rydgberg's 
series for a perfect metal and with experimental values measured on graphene layers grown on Ir 
and Ru surfaces. For free standing films, the appearance of states with binding energies in between 
the classical series is discussed. 
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Ultra thin stacks of graphene layers display a number 
of interesting properties and potential applications ow- 
ing to the lin ear bands found near the K point in the 
Brillouin zoneP^ However, to extract its full potential, 
other regions in the Brillouin Zone need to be consid- 
ered; in particular, unoccupied states in the vicinity of T 
have also attracted attention due to their potential role 
in transport of currents and heat. Indeed, the dielectric 
response of very few layers of graphene is a key phys- 
ical element to design devices based in graphene. The 
dielectric response is directly related, and therefore can 
be investigated, by looking at image states originated on 
the trapping of external electrons in the region of unoc- 
cupied states between the Fermi le vel and the vacuum 
levelP Therefore, the experimental)^ and theoreticaP 
study of image states constitute an ideal probe to better 
understand the properties of ultra-thin graphene layers. 

The image force is a non-local effect asymptotically 
dominated by correlation effectsP In order to study the 
infinite Rydberg series arising from the image potential 
one needs to compute an effective one-dimensional poten- 
tial, V(z), representing the real part of the quasi-static 
self-energy for an external unit probe charge. This self- 
induced potential is a continuos function spanning from 
inside the material, where it represents the exchange 
and correlation energy, to the vacuum region, where it 
should have the correct hydrogenic-like asymptotic be- 
havior, — i. Such a goal can only be obtained from a 
non-local spatial formalism, since a local approach results 
in a correlation potential decaying exponentially in the 
vacuum region, following the density behavior outside the 
solidP For a self-consistent first-principles theory such a 
non-local functional dependence can only be included by 
means of costly numerical calculations.- Therefore, it is 
useful and natural to search for simpler ways to obtain 
such an effective potential, which is the basic ingredient 
needed to understand the physics of image states bound 
by an ultra-thin polarizable layer like a few-layers stack 
of graphene. The simplest of these alternatives is to in- 
troduce a set of fitting parameters to continuously join 
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FIG. 1. (Color online). Self-energy, V(z), of a unit test charge 
at different positions outside and inside a thin slab (&ft = 1 
and d = 1). Black continuos line: Eq. (1). Black dots: 
RPA. Black dashed line: asymptotic law with image plane at 
zq — — — (j z |> d). Red dotted line: quartic approxima- 
tion (| z |< d) . Long-dashed and dotted-dashed (blue and 
green): Eq. (1) plus a repulsive term, Rb(z) = e~ 20(z ~ b \ Hor- 
izontal dashed and dotted lines show the first five eigenvalues 
(red and blue for even and odd, respectively); the horizontal 
thick line gives the approximate value for the work function 
in graphene. 



solutions valid either inside or outside the solid. This 
point of view has been taken, e.g., by Silkin et al. to 
study image states in free-standing graphene, joining a 
function with the correct asymptotic behavior to a local- 
density approximation (LDA) calculation for the poten- 
tial inside an atom-thick graphene layerP This approach 
makes possible to fit experimental data allowing its phys- 
ical interpretation.^ Its main weakness is its dependence 
on a few adjustable parameters, e.g. the choosing of a 
matching point in the vicinity of the surface, which can 
influence results due to the fast rate of change of the 
classical image potential near its divergence at the origin 
(image plane), etc. 
In this paper we analyze an alternative that is almost 
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parameter-free and makes a simple, flexible and accu- 
rate basis for interpreting experimental results. We use 
well-known models for the reflection of electromagnetic 
waves at the surface (infinite barrier specular modelP), 
and for the non-local static dielectric response (Fermi- 
Thomas and Random Phase ApproximatiorP) so the 
desired self-energy can be obtainedP In this approach 
only two free parameters are neeed, i.e., the electronic 
medium polarizability which is determined by the elec- 
tron density of the material, and a geometrical one give 
by the layer thickness. This approach leads in a natural 
way to a potential with proper physical features: it is 
continuous and finite over the full spatial domain and it 
has the right asymptotic behavior towards the vacuum 
region. Moreover, as for image states we are interested 
in regions in reciprocal space with k near T and ener- 
gies between the vacuum level and the Fermi energy, it is 
well justified to model graphene as a polarizable electron 
gas with a quadratic energy dispersion, as seen from the 
relevant bands for graphene and graphite for these condi- 
tions. The static dielectric response, e(k) has been mod- 
eled by a Random Phase Approximation (RPA), and by 
its small k expansion, the Fermi-Thomas Approximation 
(FT) P While the RPA yields a more accurate descrip- 
tion of excitations in the material, introducing the FT 
allows to write the potentials as analytical expressions or 
quasi-analytical ones which merely depend on a final nu- 
merical step involving the simple integration of a function 
decaying quickly for large values of the argument. The 
single parameter in these static models is the screening 
constant, kpx, that it is related to the density of states 
at the Fermi level, This value fixes the scale for 

energies, and its associated wavelength, Xft — j^;, the 
scale for lengths (atomic units are used throughout the 
paper, except where explicitly it is said otherwise). Tak- 
ing graphite as a model (2 g/cm 3 , 2s 2 2p 2 ), typical values 
for graphene are kpT ~ 1 ( r s ~ 2.5), although its precise 
value may depend on factors like doping, external poten- 
tials, etc; this is accommodated in our results through the 
scaling with kpx- The other parameter needed to char- 
acterize a thin slab is its width, 2d. For a single atom 



thick layer of graphene a reasonable value for d, should 
be related to the spatial extension of n carbon orbitals, 
dm 1. 

Self-induced potential by an ultra-thin slab. For an ex- 
ternal probe charge near a slab (Q = 1) we seek the po- 
tential acting on Q by the polarization charges induced 
in the medium by Q itself. This is obtained by comput- 
ing the total potential, and subtracting the charge's own 
naked potential. To ensure the proper boundary condi- 
tions, and according to the the specular reflection model 
at the surface, auxiliary pseudo-media are introduced for 
the polarizable slab and the vacuum that reduce the cal- 
culation to matching solutions obtained in different re- 
gions of space for homogeneous media everywhere.^ De- 
tails for the thin slab, a vacuum gap between graphene 
and a metal, and the metallic surface itself, along with 
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FIG. 2. (Color online) Bound states energies, — E n (eV), for 
the hydrogen-like series supported by the potential of the 
ultra-thin slab in Fig. [T] (rectangles, blue). These are com- 
pared with Whittaker's series (circles, red), experimental val- 
ues for Gr/Ru (triangles, purple)P and Gr/h^ (diamonds, 
green). The abscissas n labels are chosen to merely follow the 
sequence of increasing energies. 



expressions, will be given in a forthcoming publication. 
Within the FT approximation, this potential can be writ- 
ten as an expression that only depends on a numerical 
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In Fig. [T] we show the potential for a slab occupying the region — d < z < d; both in the FT approximation (black 
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FIG. 3. First four eigenfunctions for the potential displayed 
as a black continuous line in Fig.[l] Eigenvalues are in eV and 
referred to the vacuum level. For comparison, Whittaker's 
wave-functions (1,2), and the fitted harmonic oscillator wave 
function (0) are displayed (dashed). 



continuous line), and in the RPA one (black dots). In 
the all important region determining the Rydberg series 
(| z |> d), both approaches yield similar values and agree 
with the correct asymptotic power-law. Near the center 
of the slab, FT overestimates the interaction over RPA 
by about 30-40%, 
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a difference that reflects mainly in the first few states 
with a significative weight at the center of the slab. In 
units of ItFT, we see that for d > tr— , Vft(<1) = —1/3, 
and V FT (0) = -1/2 (the Coulomb holeP). Outside the 
slab, for z > Xft, the potential is well approximated by 
a classical law, 4( - z l 2o ) , corrected by an image plane, zq. 
The value of zq can be obtained by expanding the inte- 
grals for n sa 0, giving the position of the image plane in 
this model: zq — (dashed line). The potential in- 

side the slab is in turn well approximated by a quadratic 
or quartic fit (dashed red line Fig. [I]) . 

Eigenvalues and Eigenfunctions. We solve numerically 
the Schrodinger equatiorffSEHto compute the eigenvalues 
and eigenfunctions corresponding to the Fermi-Thomas 
model potential described above. Quite generally, states 
can be labelled by the number of nodes (n) , with energies 
increasing as the number of nodes increases. Further- 
more, since V(z) in Eq. 1 is symmetric, eigenfunctions 
have either even or odd parity, for even or odd n. We 
show in Fig. [I] the first five eigenvalues for V(z): dashed 
and dotted horizontal lines for even and odd parities. It 
is worth noticing the structure of this series: there is an 
isolated eigenvalue (the lowest one), while the remaining 
states cluster near the vacuum level. For standard den- 
sities (e.g., 1 < r s < 10) this first level appears below 
< —5 eV; an estimate for the work function in graphene 
(thick line). Therefore, this n = eigenstate (E = —7.2 
eV) does not fit in the standard definition for a Ryd- 
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TABLE I. Binding energies, E n (eV), compared for differ- 
ent cases. Labels n, m denote the number of nodes in wave- 
functions. Data from other authors has been labelled as in 
the original papers. We highlight in bold face numbers that 
can be compared across different calculations or experiments 
and have been given an accompanying interpretation in the 
text. 



berg state, necessarily located between the vacuum and 
Fermi levels to be observable in a standard experiment. 
Furthermore, while wave-functions for Rydberg's states 
are spatially located mostly in the vacuum region, ipo 
is symmetric and peaks at the origin (Fig. |3j continuos 
line in upper- left panel), being alike to the ground state 
of a harmonic oscillator fitting the bottom of the well 
(0.13z 2 a.u., E' = —6.2 eV), but not to states in the 
one-dimensional hydrogen-like series for a semi-infinite 
metal, that go to zero at the image plane. 

In Table I and Fig. [2] we compare eigenenergies calcu- 
lated for the FT model with experimental values reported 
for Gr/RuJ^Gr/IrPand with the limiting case of the Ry- 
dberg series for a perfect metal E m+ i — — 32 ( m 1 + i)2 ; m = 
0,1,2,.... where m refers to the number of nodes for 
each state (notice that these wave-functions only extend 
to z > half-space, and that the zero at the origin is 
not counted as a node since it derives from the boundary 
conditions). These eigenvalues can be conveniently ob- 
tained from multiple-scattering techniques^ while wave 
functions are obtained as the solution to Schrodinger's 
equation with Vjj (z) — — (z > 0), that can be reduced 
to Whittaker's differential equation.^ The similarity of 
values found for the antisymmetric (n~ ) members of the 
series of states for V(z) and Rydberg's series is striking. 
Such a similarity can be better understood by looking 
at the corresponding wave- functions (Fig. [3]). Boundary 
conditions make all Whittaker's wave-functions to go to 
zero at the origin, a condition that in the case of a sym- 
metric well can only be fulfilled by odd wave functions. 
Moreover, if n~ and m give the number of nodes for odd 
wave-functions for the symmetric potential, and the Ry- 
dberg's one respectively, we can make a one-to-one cor- 
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Whittaker's 12 3 4 

z 3 12 29 51 79 

Eq. (1) 1 2 3 4 5 6 7 8 9 

z 3 6 12 18 28 34 48 58 78 



TABLE II. Expectation value z (A) for Whittaker's and 
Eq. (1) wavefunctions. 



respondence, - 9 ~ 1 = m, that simply tells us that both 
sets of wave-functions have the same number of nodes if 
wT is divided by two (only half-space) and the node at 
the origin is discounted. From a physical point of view, 
we can envisage two relevant limits: a free standing slab 
producing a symmetric potential with states labeled by 
the number of nodes and their parity, and a slab lying on 
a substrate where a particular surface gap may prevent 
penetration of wave- functions inside the material leaving 
only half the space accessible for image states. Therefore, 
it is reasonable to argue that, according to our analysis, 
for a free standing slab, or one interacting weakly with 
the substrate support, it should be possible to observe the 
even states as new energies located in between the usual 
hydrogenic ones. It is interesting to notice that one of 
these states may have been observed for Gr/Ru (n = 2,)p 
while none of these have been reported for Gr/IrW This 
fact must be related to the strength of the interaction 
between the supporting metal and the graphene layer, 
and the penetration of wave-functions in both metals in 
a way that goes beyond the scope of the current analysis. 

To assess how sensitive are the eigenvalues to the de- 
tails of the model potential we have added to V(z) a 
repulsive term modeled as an exponential wall: Rb(z) — 
e -20(z-6)_ Q n a me t a lli c surface such a "repulsive" term 
can originate because electronic gaps existing for partic- 
ular surface orientations. The resulting potential for the 
repulsive barrier located near the slab surface (6=1, 
green dashed line in Fig. 1) is similar to the classical se- 
ries with an image plane to prevent the divergence at the 
origin ( 4 ( z + Zo ) ) arid can be understood by introducing a 
quantum defect in Rydberg's series. The barrier, on the 
other hand, can be introduced below the surface (e.g., 
b = 0, blue dashed-dotted line in Fig. 1), and a numeri- 
cal solution for Schrodinger's equation can be obtained. 



Results in Table I show that, for large m, eigenvalues 
tend to the classical ones in accordance with the fact that 
the average position of the electron is far away from the 
surface, therefore mostly influenced by the asymptotic 
behavior. Table II gives the expectation mean values in 
A, z = J °° ip(z) zi[>(z) dz, for Whittaker's wavefunctions 
compared with the ones corresponding to Eq. (1). These 
values compare well, which reflect the manifest similar- 
ity between wave functions commented on Fig. [3] The 
fact that z > d for u > 1 determines a spatial region 
where both FT and RPA dielectric functions lead to ap- 
proximately the same potential, suggesting that higher 
k-corrections to the dielectric function arising from the 
random phase approximation are negligible, at least for 
n ^> 1 states. Taking away the first level, largely af- 
fected by the details near the bottom of the potential, 
the rest of the series is only modified by a percentage 
comparable to differences found in Table I between sim- 
ilar entries. Therefore, we argue that our assignation of 
levels is sound and robust. 

Conclusions. Using standard models for the dielectric 
response and the reflection of electromagnetic waves at 
a surface we have computed the static self-energy for an 
ultra-thin slab mimicking a graphene layer. The self- 
induced potential goes continuously from the exchange 
and correlation energy inside the material to the classi- 
cal asymptotic image potential in the vacuum. Eigenval- 
ues and eigenf unctions have been compared with Whit- 
taker's classical series and recent experiments on Gr/Ir 
and Gr/Ru. The odd members of the series for the 
slab show a remarkable resemblance to the solution of 
Schrodinger's equation for the classical image potential 
(Whittaker's wavefunctions). On the other hand, even 
wavefunctions arise as new states that differ from Whit- 
taker's in several key respects, e.g. their non-zero density 
probability at the origin. For the case of films weakly 
interacting with a support some new states may conse- 
quently appear in between the classical ones, that can 
be traced back to the even states in a free-standing slab. 
Such a case could have been observed in recently mea- 
sured experimental values on Gr/Ru. 
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